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Abstract: The Horndeski theory is known as the most general scalar-tensor theory with second-order 
field equations. In this paper, we explore the bi-scalar extension of the Horndeski theory. Following 
Horndeski’s approach, we determine all the possible terms appearing in the second-order field equations of 
the bi-scalar-tensor theory. We compare the field equations with those of the generalized multi-Galileons, 
and confirm that our theory contains new terms that are not included in the latter theory. We also discuss 
the construction of the Lagrangian leading to our most general field equations. 
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1 Introduction 

The inflationary scenario [1] is now regarded as a necessary ingredient of modern cosmology because 
the recent observations of cosmic microwave background anisotropies [2-5] strongly suggest the epoch of 
inflationary expansion in the early Universe. Since inflation is typically driven by scalar field(s), scalar- 
tensor theories provide a firm framework to study the dynamics of inflation. On the other hand, the 
unknown energy called dark energy is shown to be dominant in the present Universe [6-8] and might be 
understood as an outcome of infra-red modification of gravity. Scalar-tensor theories provide powerful 
tools to realize such a modification, and their phenomenological features have been studied intensively 
to confront them against observations. Although a plethora of inflationary models and modified gravity 
theories have been proposed thus far, unfortunately, we have not yet succeeded in finding the real theory 
and have kept seeking for it. We have two options to address this problem. The first one is to pursue the 
ultimate (real) theory on the basis of theoretical consistencies, which is often called top-down approach. 
The other is to construct a framework of theories as general as possible, which is called bottom-up 
approach. Of course, though both of approaches are complementary, one of the merits to take the latter 
approach is to give a unified understanding of various models proposed individually. Another is that one 
can easily pin down or constrain models once characteristic observational results would be reported. 
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The Horndeski theory [9] provides a typical working example of the latter approach because it is 
the most general single scalar-tensor theory with second-order field equations. It is shown [10] that the 
Horndeski theory is equivalent to the generalized Galileon [11]. In fact, almost all of the inflationary 
models with single inflaton proposed so far can be described by this theory in a unified manner. Various 
aspects of single field inflationary models have been studied in this framework [10], which is useful for us 
to constrain the models from observational results. 

In this paper, we take the latter approach and try to extend the Horndeski theory, which includes 
only one scalar degree of freedom, to the bi-scalar case.^ Although no observational results suggest the 
presence of multiple (light) scalars during inflation, our framework can clarify the essential difference 
between single and multiple scalar cases, and is useful for further constraining or detecting the presence 
of multiple scalar degrees of freedom from observations. Such an attempt to extend the Horndeski theory 
has already been discussed. First of all, multi-field Galileon theory was proposed in the flat spacetime 
[17-20]. Later, the covariantization of this multi-field Galileon theory, called generalized multi-Galileon, 
was considered [21] and conjectured that the theory would correspond to the multi-field extension of the 
Horndeski theory, that is, the most general multi-field scalar-tensor theory with second order equations 
of motion. Though the multi-field Galileon theory in the flat spacetime is proven to be the most general 
multiple-scalar field theory in the flat space-time with second order scalar equations of motion [22] , it was 
shown that the generalized multi-Galileon is not the most general theory [23] because this theory does not 
contain the multi-DBI inflation models [24-28] , in particular the double-dual Riemann term appearing in 
these models. Motivated by these considerations, in this paper, we try to construct a multi-field extension 
of the Horndeski theory. Especially we focus on the extension to the bi-scalar case as a first step. 

The organization of this paper is as follows. In Sec. 2, we construct the field equations of the most 
general two-scalar tensor theory with second order field equations following Horndeski’s procedure. This 
section is the main part of this paper. In Sec. 3, we compare our theory with the generalized multi-Galileon 
theory. We show that terms which are missing in generalized multi-Galileon theory are actually contained 
in our theory. In Sec. 4, we comment on the construction of the Lagrangian corresponding to the held 
equations we obtain. Finally we summarize our paper and discuss the results in Sec. 5. 


Notations and conventions 


Before closing the introduction, we summarize the notations and conventions used throughout this paper. 
We consider a four-dimensional spacetime with a metric gab and two scalar helds (p^ with I = 1,2. 
Following Ref. [9], derivatives of gab and with respect to the coordinates are denoted as 

I _ dp^ 


9ab,c — 


d^ab 
dx^ ' 


dx^ ’ 


( 1 . 1 ) 


respectively. We denote the covariant derivative of p^ with respect to gab and its scalar product respectively 
as 

Pi ^ VaP^, ( 1 - 2 ) 

^Another direction is to consider a scalar-vector theory instead of a scalar-tensor theory. Motivated by the earlier work 
of Horndeski [12], construction of the most general vector theory with second order field equations on flat space, called the 
vector-Galileon theory, was attempted in Ref. [13] . Recently, the scalar-tensor theory with higher order equations of motion 
without introducing the ghost was proposed as well [14, 15]. See also ]16] for yet another possible way of extending the 
Horndeski theory. 
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where is symmetric in I and J. We use a strike “ | ” also as a separator in (anti-)symmetrization. 
For example, [I\JK,L\M\ stands for anti-symmetrization of I and M. Partial derivatives of a function 
dg, d'^g, are expressed as 


j^a...b;cd _ 


dA 


a...b 


dgcd 

dA^-b 


j^a...b\cd^e _ 


dA 


a...b 


Aa-b; ^ 


dgcd,e 

dA^-b 


Aa...b\c _ 

d(t>{ 


j^a...b\cd^ef _ 


dA 


a...b 


A 


dgcd,ef 

Qj^a...b 


a...b]cd _ 


(1.3) 


and partial derivatives of a function A{(j)^, X'^^) are expressed as 

, _ dA _1 f dA \ 

“ d(/)^' “ 2 dX'^^) 


(1.4) 


In the equations of motion and the Lagrangian, we use the generalized Kronecker delta defined by 




xIK _ nr/ s:K 
OjL = . 


Repeated indices are summed over o = 0,1, 2, 3 and I = 1,2. 


(1.5) 


2 Construction of the most general equations of motion 

The first step of the construction of the most general scalar-tensor theory of Ref. [9] is to work out the 
most general equations of motion that are of second order in derivatives and compatible with the general 
covariance. In this section, we generalize this construction to the case with two scalar fields. 

2.1 Assumptions 

The assumptions imposed on the theory we are going to construct are summarized as follows. 

1. The theory has a Lagrangian scalar density, C. 

2. The Lagrangian scalar density, C, is composed of a metric, two scalar fields, and their derivatives 
up to arbitrary order: 

^ ^ (^9abj 9ab,cj 9ab,cdj ■ ■ ■ j 4* j 4^,ab^ ■ ■ ■ ) ; (^-1) 

where / = 1,2. 

3. Field equations are composed of the metric, the two scalar fields, and their derivatives up to second 
order: 

bH, 

0 ~ T ~ ^/~9 G {dcd’i 9cd,ei 9cd,ef': 5 4*,cd) ? (2*^) 

^9ab 

6/1 

0 ~ ~ ~9 {9abj 9ab,ci 9ab,cd] 4^ i 4^,aT 4^,ab') ) (^-3) 

where 5Cj6A is the variation of £ with respect to a field A. 
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Let us consider the variation of the action under an infinitesimal coordinate transformation, ^ 
, which is given by 

6 jd^xL = 2 j (2.4) 

It follows from the assumption of C being a scalar density that Eq. (2.4) vanishes identically, implying 
that the integrand itself vanishes since may be chosen arbitrarily. Thus, the identity 

(2.5) 

holds. In the case of pure Einstein gravity this identity reduces to the well-known contracted Bianchi 
identity, = 0. In this sense, the identity (2.5) may be regarded as a generalization of the Bianchi 

identity. We are trying to construct the most general bi-scalar-tensor theory with second-order field equa¬ 
tions using the identity (2.5). Because both G°'^ and £j are assumed to be of second order in derivatives, 
the left-hand side of Eq. (2.5) would yield third derivatives in general while the right-hand side contains 
at most second derivatives. This indicates that the left-hand side of Eq. (2.5) must be free from third 
derivatives. 

The construction of the most general second-order equations of motion for bi-scalar-tensor theories is 
divided into two parts: we first determine the most general second-order rank-2 tensor whose divergence 
remains of second order. After that, we impose the identity (2.5) on the rank-2 tensor to constrain its 
form. 


2.2 Second-order rank-2 tensor whose divergence is of second order 

In this subsection we construct the most general second-order rank-2 tensor, 0“^, whose divergence is also 
of second order. The conditions that S/aG°'^ has no third derivatives can be expressed as 


d(jcd,efg 

(2.6) 

■ 

(2.7) 

Using the chain rule, VbG°'^ is rewritten as 


V 7 r>ab f->ab:cd,ef ^ , ^ab:cd,e^ , f->ab\cd^ , ^ab;cdil . ^ab-,c il . ^ab; si , ^bc-r^a , ^ab-r^c 

VbG =G ’ ' gcd,efb + G ’ ' 9cd,eb + G gcd,b + G j cp^cdb + y i4>,cb + y i9,b + y ^bc + y Uo 


(2.8) 

where T^^ are the Christoffel symbols. With the help of Eq. (2.8), we can 
and (2.7) are equivalent to 

show that the conditions (2.6) 

Qab;cd,ef _j_ Qae;cd^fb _|_ Qaf;cd,be _ q 

(2.9) 

pfabicd 1 pyacidb , pyadihc n 

y I +y 7+^ 7=*J’ 

(2.10) 

respectively. From the “invariance identity” (see Refs. [29, 30]), we have 


gab;cd,ef _|_ gab-,ce,fd _|_ gab-,cf,de _ 

(2.11) 
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( 2 . 12 ) 

(2.13) 


By repeated use of Eqs. (2.9), (2.10), and (2.11), we obtain 


Qab-,cd,ef _ Qcd\ah^ef _ Qef;cd. 


ab 


r»<ib\cd r>c.d:ab 

Q I =G j . 


For convenience, we now introduce the notion of property S following Ref. [9] . A quantity ^“i“ 2 ---a 2 n-ia 2 n 
is said to have property S if it satisfies the following conditions: (i) it is symmetric in {a 2 e-i,a 2 e) for 
i = 1,2,... ,n; (ii) it is symmetric under the interchange of any two pairs { 021 - 1 , a 2 i) and {a 2 m-i, 02 m) 
for £,m = l,2,...,n; (hi) it vanishes if any three of four indices, { 021 - 1 , 021 ) and {o 2 m-i, 02 m) for 
i,m = 1,2,... ,n, are symmetrized. It is shown in Corollary 2.1 of Ref. [9] that ^“i“ 2 ...a 2 Ti-ia 2 n van¬ 
ishes if ^“i“ 2 ...a 2 n-ia 2 n property S and n > 4 in four-dimensional spacetime. 

Let us introduce the quantity i?“i“2...a4n+2m+ia4n+2m+2 defined by 


2^aia2...a4n+2m + ia4n + 2m+2 = jjy ^ j 


d 


dga^i- 


a4i-i an ,a4i+ia4i+2 


|n™rl 


d 


■>aia2 




“4n+2j + l<l4Ti+2j+2 . 


(2.14) 


This is a n-th derivative with respect to gaiai+ 4 ,ai+ 2 ai +3 m-th derivative with respect to 4>^^aiai+i 
Using Eqs. (2.9)-(2.13), one can easily check that ^“i“2...a4n+2m+ia4n+2m+2 j^g^g property S. Then, Corollary 
2.1 of Ref. [9] implies that B°-i<^2-a4,n+2m+ia4n+2m+2 vanishes for 2n-|-m > 4, leading to the following three 
sets of identities: 


dgcd,ef dgij^kl 

d d _ 

dgcd,ef dcj)^- 
d d d d 


(2.15) 

(2.16) 
(2.17) 


By integrating Eqs. (2.15)-(2.17), we can determine the form of the gravitational field equations. 
First, integrating Eq. (2.15) yields 


gab ^ ^ ^ab ^ ^abcdef ^ ^ab ^ 


(2.18) 


where and are functions of gab,gab,c,4>^,4>^^a^ 4’^ab- Note that we have used the identity 

e^'^^^gcd,ef = Recdf + (2.19) 

at the second equality of Eq. (2.18), where are functions of gab, gab,c,4>^,4'^a It can be seen 

that and have property S. Substituting Eq. (2.18) into Eq. (2.16) and integrating it, we obtain 

gab ^ ^ Rcdef + (2-20) 

where ^abcdefgh functions of gab, gab,c, 4>^, and cj)^a^ while are functions of gab, gab,ci 

(ji^a and (j)^ab- Here again it can be seen that ^abcdefgh ^ab j^gve property 5. Repeating the 

same procedure and integrating Eq. (2.17) give 

e“‘ = (2.21) 


- 5 - 







where all of the above ^ tensors are composed of gab-igab,c-,(t^^■, and and have property S. Although 
our final goal is to determine the most general equations of motion for the bi-scalar-tensor theory, the 
equations given up to this point hold irrespective of the number of the scalar fields. 

Our remaining task in this subsection is to construct explicitly all the possible ^ tensors that have 
property S and are composed of gab,gab,c,4’^, and For this purpose, we can use g^^, and the 
totally antisymmetric tensor as building blocks, from which the ^ tensors are built by taking their 
products and linear combinations appropriately. There is no elegant way, and what we will do is to 
exhaust all the possible combinations of those building blocks yielding the ^ tensors. Let us begin with 
the simplest one, It is not difficult to find that the following one is the most general symmetric rank-2 
tensor composed of , 5 “^, and g(“Lc. 

\a 

( 2 . 22 ) 


where a{(j)^, X'^^) and {(p^ X"^^) are arbitrary functions of (p^ and X"^^, and has the sym¬ 
metric property 6 ^*^ = Here, we have used, for the first time in this derivation, the assump¬ 

tion that the number of the scalar fields is two, which greatly simplifies the expression of ^ tensors 
and the following procedure. Without this restriction we would have for example the term such as 
cijKL ™ where cijkl is arbitrary functions of cp^ and X^-^, 

and totally anti-symmetric in /, J and K. In a similar manner, we work out 


^ (^gacgbd ^ gadgbc _ 2g<^bgCd^ 

+ bijK - 2 

+ CIJKLM 

+ dijKLM , 

where a/, bijK, cjjklm, and djjKLM are arbitrary functions of (p^ and X"^^ satisfying 


(2.23) 


bijK = biKj, 

CIJKLM = CjkJLM = CijkML = CjlMJK, 

diJKLM = —diKJLM = —diJKML = duMJK- (2.24) 

The explicit forms of ^abcdefgh ^ ^abcdefgh bi-scalar case are given in Appendix A. 

Substituting all the ^ tensors into Eq. (2.21) and rearranging the equation, we arrive at the most general 
second-order rank -2 tensor 0°“^ whose divergence is also of second order, 

dt = ASt + + C,S»‘!f + D,jKSSg<PlJ-''Ul^'' + E,jKLM<>^<i>[yui!<P‘-''</'‘p'' 

+ EuKLuiZil 

+ + liZSRj + 

+ (2.25) 

^One might consider other rank-2 tensors such as but this tensor is excluded because it is not symmetric in 

a, b. 
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where A, Bij,Ci, Dijk, Eijklm, Fuklm,Gij, Hjjkl, I, Jij, Kj, and Ljjk are arbitrary functions of 
and , and they are subject to 


Bij = Bjj, 

Dijk = 

Gij = Gji, 

HijKL = 

II 

Lijk = 


Djik, 

= Hjikl = Hulk, 
Ljik = Likj- 


EijkLM = —EkJILM = —EjlKJM = EjjlkM, 
EijkLM = —EjikLM = —EijLKM = FjilkM, 

( 2 . 26 ) 


2.3 Consequence of Eq. (2.5) 

In the previous subsection we have obtained the most general second-order rank-2 tensor whose divergence 
remains of second order, Q°‘^. As can be seen, involves many arbitrary functions. It turns out, however, 
that those functions are not completely independent in order for the equations of motion to be compatible 
with general covariance. In this subsection, we impose the identity (2.5) that arises due to general 
covariance, i.e., we require that the divergence of 0°“^ is written as a product of 0^1“ and some scalar 
function. This procedure will reduce the number of the arbitrary functions. A straightforward calculation 
shows that 

+ + f^lS^bdRcl “</>''' + 

+ 2FUKLM + 

+ ^^bdfh (^FijkLm 4>^^°'4>'^^^'^ (^FijKLM4’\c4'\e4'^4'fh) *^1^'^’ (2.27) 

where the coefficients are functions of (j)^ and and defined as 

a/./ =Gu - 2 Jij + 2 Ki^j - 2 HklijX^^, 

(3ij = - I,ij + Jij - Kjj + 2 Jkl,ijX^^ , 

lIJKL = — 2Jij^kl + Hjkjl, 

3 

eijK =K{i,k)j - 2^kij, 

7^/ ~2^i + ‘^^,1 ~ {Hjki + 8Jj[k,i]) + AEjklmiX'^^X^^, 

I'l.jKL = - Gij,kl + ^Hk(ijl) + 2Hmnij,klX^^ - 3Luj,k, 

^IJK = - Gi^jk + 2Dj^ik) - 2Gik,j + 2 {Dlmi,jk - 4Hy\uK,\M]) X^^ 

1 ^ T71 ^LM Q TP \^LM 'y'NO 

- - otjLMNOI.JK^ ^ , 

iiJ = - Xjj + Bu - Gjj - 4 .Dk[i\j^\l]X^^ 

Q jp \rKL 'W'MN I 1 ^ IT' 'yKL \rMN 

- OI^kLMNJJ^ + ^^^KIMNJ.L^ ^ , 

CijK = - 2^UK - 2 Jij,k + 4ElmijkX^^, 
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VIJKL =-^HjjkL, 

^ijKLM =-Dijk,lm + Hijkm,l “ ‘^Emlijk “ ^Enoijk,lmX^'^, 

C^IJKLMN =HijkL,MN “ HimNL,JK, 

TIJKLM = - L[jijk,L\M], 

I'lJK = - (2.28) 


We present the explicit form of Qj in Appendix B, though it is irrelevant to the following derivations. 

In order for the right-hand side of Eq. (2.27) to be proportional to all the terms that are not 
parallel to must vanish identically. We now derive the conditions for this following the procedure 
of Ref. [9]. Let us first focus on the ejjx and ljjk terms in Eq. (2.27), which are proportional to 

,/|a ,J|c„ eg. 

^\b ^\d ^fh ■ 


The coefficient of 4*^ ,op9qr,st iri this quantity can be extracted by taking a derivative with respect 

to (p ^rnn4 ,op9qr,st 




+ 2iijK5Z;i 6^^9^'^^g^'^f. 


Equation (2.30) must vanish when contracted with a vector Ya such that = 0 because Eq. (2.5) 

implies that 




(2.31) 


Thus, we obtain a constraint equation given by 


26(a| 

+ 2iijKYJZfl 6^/g^Y = q. (2.32) 

This constraint has eight free indices m, n, o,p, q, r, s, t, and its any component must be fulfilled. We first 
take the trace of Eq. (2.32) by contracting with gmngop9qr, giving 


- 8eyA\j\B)<P^^^Y^^ + 8qA\J\B)4>^^^Y^^ = -8eyA\J\B)(^^^^Y^^ = 0 , 


(2.33) 


^ To obtain this result, we use the fact that derivatives of (I)i,ab4>j,cd and Rated are given by 

d \abP |c<i) j./ r-m, r-n, r-o cv , cJ cm c-n c/ ro c-p ^ / r/ cJ j-\mnop , cJ cl j-\mnop\ 


dRabcd 1 / ■p.qrst , j-\qrst j-\qrst 

Qg ^ ^ \^abcd ~r ^edab ^abdc ^bacdj ’ 


where 
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where we have used i[A\j\B) = 0 which follows from the definition (2.28). This equation must be satisfied 
for any and therefore it is necessary to impose e(^A\j\B) (= ^ajb) = 0. Further constraints can be 

derived from Eq. (2.32) as follows. Let us project Eq. (2.32) to the basis vectors (7 = 1,2) 

that satisfy 

YaY^ = YaY^ = 1, Ta^“ = 0, = !> = 0 for 7 = 1,2. (2.34) 

By contracting Eq. (2.32) with where and are either y“ or y“, we 

find 


0 = -Ae^A\j\B)^bdYay^{ZcW'^ + WcZ'^) 

+ AiAJB5MZaV''{ZcW'^ + WcZ^)X^^ (^x^D\{FJ^G)\E) 




(2.35) 


Using e(^A\j\B) = 0 obtained in the previous step, we find lajb = 0 as another constraint to be imposed. 

Repeating a similar procedure for any other products of the second derivative terms in Eq. (2.27), we 
find the following constraint equations for the coefficient functions: 

= -2/3j/, '^IJKL = -4:'ni[J\L\K], OtAI ^I{cZ^D)B - ^IbXcD = ^ (^(CD)ylS “ ^(C|SA|D)) j (2.36) 
(^IJK = I'lJK = ^IJK = \lJ)KLM - \l\LM\J)K = fJ-I = Cl[JK] = = TUKLM = FjjkLM = 0, (2.37) 

FACDbX^P — Fab(E\cX\f)D ~ ‘^<^{EF)CADB = 0; ^B[A\K\C] = (^EFICA]DB = 0, (2.38) 


where XjJ is the inverse matrix of X^'^. 

The constraints (2.36) and (2.37) impose the following conditions on the functions appearing in 
Eq. (2.25): 


Bu = -2 (y + 2>V)^,,J + Ajj + - 16Ek^hmn\j),lX^^X^^ 

- 8 {Jk(i,j),l - Jkl,i,j) X^^, (2.39) 

Cl = -2 (y + 2>V)_, + 2 {Djki + x-^^ - SEjklmiX'^^ X^^, (2.40) 

Fijklm = 0, (2-41) 

Gij = 2Jjj - 277(,,j) + (2.42) 

HijkL = ‘^JiJ,KL: (2.43) 

%.j] = -2Jk[i,j]lX^\ (2.44) 

Ki,jk = Kjjk, (2.45) 

Lijk = (2.46) 

7=iy + W, (2.47) 

where >V = >V(<^^), and F = ,X^^) is a function satisfying Fjj = G/j, which is integrated to give 

F = j Gij dX^J = I {2Jij - 277,.j + AJkijlX^^) dX^-^. (2.48) 
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The conditions C,i[jk] = ^IJK = \ij)klm - \i\lm\j)k = 0 in Eq. (2.37) imply 

Di[jk] = -‘^Ji[j,K] + ■, (2.49) 

I^I(JK) = 2^J,IK + Gjk,I + {—DlmJ,IK + 1 ^]) 

+ 8E^jiilm\k)X^^ + 4.Elmnoj,ikX^^X^o, (2.50) 

0 = 2^IJK,LM - 2^{I\LM,\J)K + HijkM,L - -f^(7|LMA',|J) 

— 2 [Eml{ij)k — Ek{ij)lm) — 4 {Ejyfo{ij)K,LM — Eno{i\lm,\j)k) , (2.51) 

and Eq. (2.38) implies 


Giij,k]l = 0, (2.52) 

Hijk[l,m]n = Oj (2.53) 

G(ij,kl) = + 2Him{ij,kn)X^^ - 2K(^i^jk)^l ■ (2-54) 

Equation (2.52) is nothing but the integrability condition which guarantees Ejj^kl = X^kljj, and hence 
the integral (2.48) indeed exists. 

2.4 The most general second-order equation of motion 

Now we are at the final stage of deriving the most general second-order field equations of the bi-scalar- 
tensor theory. Substituting Eqs. (2.39)-(2.47) into Eq. (2.25), we at last obtain 

Gt = A6^ + [-2Ej - m,i + 2 {Djki + 8Jj[k,i]) X^^ - SEjklmiX^^X^^] 

+ {-2Ej,j - 4>V,/,j + Ajj + 2Djkj,lX^^ - 16Ekimnj,lX^^X^^ - 16Jk[i^luX^^) 

(2.55) 

This is the main result of this paper. The most general held equations for the single-scalar case [9] are 
reproduced as should be if one restricts the number of the scalar helds in Eq. (2.55) to one. Note that 
one can eliminate W(0^) from the above equation by redehning E —?■ E{cj)^, X'^^) = E + 2W’. We can 
see that Eqs. (2.50)-(2.54) do not reduce the number of the arbitrary functions because these are the 
relations between derivatives of the functions. In other words, these do not affect the structure of the 
held equations (2.55). As we will comment in the hnal section, these may, however, help us to check the 
integrability conditions for the held equations. 

Erom the relation (2.5), the scalar-held equations of motion are found to be 

£l = 2Qi + V|^• (2.56) 
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3 Comparison with the generalized multi-Galileon theory 


The covariant version of the multi-Galileon theory in the flat spacetime [21] was conjectured to be the 
most general multi-scalar-tensor theory with second-order field equations. However, later it was pointed 
out that this theory is not the most general one [23]. A counter-example is given by the multi-field DBI 
Galileons [25]. In this section, we compare the most general second-order field equations obtained in the 
previous section with the field equations of the generalized multi-Galileons, and identify the terms that 
are missing in the latter theory. 

The action of the generalized multi-Galileons is given by [21] 


1 

7 ^ 


£ 


G2 - + G^R + G4,/j 



(3.1) 


where G 2 ,G^i,Gi and G^i are arbitrary functions of 4>^ and , and is the Einstein tensor. The 
functions G^j^jk, G^jj^kl, G^j^jk and G^i^jk,lm are totally symmetric with respect to all of their indices, 
/, J, AT, L and M in order for the field equations to be of second order. A straightforward calculation leads 
to the held equations for the generalized multi-Galileons, 

E ^\ C ) = 5“' + + 03(7,j) - ^4,7,7) 

+ {-X^^Gsijk + G4,7 + 2X-^^G^u,k) 

+ (^--^G^uk + 2G^k{i,j) - -jG^K,i,^ 

_ 1 (G 4 - 2G47jX^-' + G5(7,j)X^'') 

+ + ^^^G^ijkl - -jG^f^i^j) - -X^^G^ijk,!^ g’'^'^^’idr^\c'^^\y 

+ ^ (G 47 J - G5(7,J)) 

+ 2 [Gaijkl - G^kl{i,j)) 

- i (G57J7. + X^^G.uklm) (3-2) 

Comparing Eq. (3.2) with Eq. (2.55), it is easy to see the exact correspondence between each term. It is 
also found that terms corresponding to Ejjklm are lacking in the generalized multi-Galileon; this is a 
completely new term. We would, however, emphasize that, even setting Ejjklm = 0, Eq. (2.55) covers 
a wider class of theories than the generalized multi-Galileons. This fact is to be illustrated in a concrete 
example presented below. Note in passing that the coefficient functions of the above equation satisfy all 
the constraints (2.36)-(2.38) found in the previous section. 

The double-dual Riemann term deduced from the multi-held DBI Galileons, 




(3.3) 
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is not included in the Lagrangian of the generalized multi-Galileon theory [23]. One can however check 
that this term is actually contained in our theory. It is straightforward to derive the field equations from 
Eq. (3.3): 

This is reproduced by setting Jjj = 2 {6ij5kl — ^ik^jl) in our field equations. 


4 Candidate Lagrangian 


Having thus determined the most general second-order field equations of the bi-scalar-tensor theory, let 
us now explore the Lagrangian that gives the field equations we have derived. For the construction of the 
Lagrangian, we employ the same strategy as taken in Ref. [9]. In the single scalar-field case, Horndeski 
found that the form of the Lagrangian can be guessed from the trace of the gravitational field equation. 
In the same way as in the single-field theory, we take the trace of the field equations of the bi-scalar-tensor 


theory (2.55) and arrive at the terms of the following form as a candidate Lagrangian: 

£i = (4.1) 

C2 = V^ , (4.2) 

= ^ (4.3) 

£4 = . (4.4) 

A = , (4.5) 

= (4.6) 

where Mj'^\ Mj^J, and arbitrary functions of 4>^ and satisfying 

M?)k = v'L. (4.8) 

Mf] = Mf], (4.9) 
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In order to maintain the second-order equations of motion for the scalar fields, we have to impose extra 
conditions on these functions. For example, the Euler-Lagrange equation of £1 for the scalar field cj)^ is 
given by 


+ - 2V7.LA'"'' + 




^J,KI^\c ^\d 


^I,JK^\c ^\d 


(4.11) 
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The two terms in the first line of Eq. (4.11) are of third order, while the other terms are of second or first 
order. To eliminate the third-order derivatives, we therefore impose 


= »■ ('‘■12) 
Performing the same analysis, we find that higher-derivative terms are removed by requiring that 




— 0 


= 0 , 




— 0 

I,J[K,L]M 


M 


(5) 


IJ,K[L,M]N 


= 0 , 


M 


(7) 


IJKL[M,N]0 


= 0 . 


(4.13) 


The Euler-Lagrange equations for the Lagrangian densities (4.1)-(4.7) are listed in Appendix C, from 
which one can see the relations between the functions appearing in Eq. (2.55) and those in Eqs. (4.1)-(4.7): 


.4 = + iMfjX'-’ - - 8 - M^j,, 






(4.14) 


Dijk = + ^^(7^) + ( Mr/r n ...r.r.-M 


r(3) 

^L{IJ),MK 


r(3) 

‘(IILM,IJ)K 


(3) 

IJL,KM 


X 


LM 


- 2M 


(4) 

K,I,J 

r(7) 


-h2 2M 


r(5) 


-M 


(5) 


IJ,K 


r(5) 

L(I,J),MK 


r('i) 


X 


LM 


-LIOA#)'-' _1_ C A/fO; Y-LjIVI y- 


(7) 


MM vNO 


(4.15) 


EijkLM — - 


— 4M^^^ -I--I- 

2 PRQ,SM ^^^^PR,Q,SM ^ ^^^PRQSM ^ ^^^PRMSQ ^ ^'^^PRQMS 


;(7) 

MRQSN,OM 


r(7) 

MRNSQ,OM 


r(7) 

MNQSO,RM 


X 


NO 


T + 2yV = - 2M^f]x^‘^ - 2Mf]x^'^ -h 2Mf)x^-^ -h jjX^-^X^^, 


M,J 


MK,JL^ 


1 


r(2) 


Ju = + M),’ + ( 2M 


(4) 


r(5) 


2 

Kl = -Mfi-X''- 


MJ 


r(5) 

M{I,J)L 


-M 


(5) 


IJ,KL 


X 


KL 


(4.16) 

(4.17) 

(4.18) 

(4.19) 


In addition, comparing the and f^rms, we see that the following two conditions must be 


satisfied: 


A 


KL 


- 2 (A + 2>V)^,^^ + Ajj + 2DKiij),LX^^ + WEkmnujuX^^X^^ - 8 {Jk{i,j),l - Jkl,i,j) X 

= M,W , + 24/^1, - + 2M^Xl.J> - 2<’(b),l) 

- 4 (vfAi + <lu - V'""' + jT 

_i_ Q ( _9 -I- 

O y^J-MNKL{I,J) ^^'^^MN(KI)J,L ^'^^MN{IJ)K,L 


ML 


r(6) 

,/J 


XKLx^N 


(4.20) 


- 2 (A + 2 W )^^ + 2 {Dijk + ^Ji[j,k]) " SEijlmkX^^X^^ 

= -Xli%X^^ - 2 (m® + 2Mi%X^-^) + + 2M®^ (4.21) 
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The Lagrangian is constructed by solving the above equations for , M^’^\ and unfortunately 

we have not accomplished this step yet. In Appendix D, we review the construction of the Lagrangian for 
the single-scalar case. The lesson from the single-field Lagrangian is that we should probably integrate 
Eq. (4.19) first to identify compatible with Eq. (4.13). We have not yet succeeded even in solving 
these equations. In addition, it should be kept in mind that the terms we considered 
might not be enough to construct the true Lagrangian. Those terms are simply inferred from the trace of 
the equations of motion, and in general there is no guarantee that they are enough though they happened 
to be so in the single scalar-field case. Albeit these difficulties, we believe that above calculations are 
useful to build the Lagrangian for our theory. We hope to report on this hnal part of the construction of 
the most general second-order bi-scalar-tensor theory in the near future. 


5 Discussions and summary 


In this paper, we have reported our attempt to construct the bi-scalar generalization of the Horndeski 
theory in four-dimensional spacetime. Eollowing Horndeski’s method, we have succeeded in deriving 
all the possible terms appearing in the most general second-order field equations for the bi-scalar tensor 
theory. We compared our held equations with those of the generalized multi-Galileon theory, and identihed 
the terms that are really not included in that theory. In particular, we conhrmed that the double-dual 
Riemann term, which was shown to be missing in that theory [23], can be reproduced from our results by 
choosing the arbitrary functions appropriately. We have also discussed the construction of the Lagrangian 
yielding the held equations we found. For the construction, we have taken Horndeski’s approach based 
on the trace of the held equations as a candidate Lagrangian. From the Euler-Lagrange equations we 
have obtained several differential equations for the functions in the Lagrangian, though we could not solve 
them to give explicit forms of the functions. In fact, it is still unclear whether or not the candidates of 
the Lagrangian we proposed suffice to generate all the terms in our most general held equations. 

As discussed e.g. in Ref. [13], we have to impose the integrability conditions on the held equations in 
order to ensure the existence of a corresponding Lagrangian. The integrability conditions are summarized 
as 



Sy/^G'^^iy) 

Sgpxiy) 

Sgpuix) 


SV^£i{y) 

S4>^{y) 

Sgpu{x) 

Sy/^£l{x) 

^y/^£j{y) 

5(j)J{y) 

5(t)\x) 


(5.1) 

(5.2) 

(5.3) 


where 5/5A denotes variation with respect to a held A. Surprisingly, Horndeski found the corresponding 
Lagrangian in a rather heuristic way without using the integrability conditions explicitly. This implies 
that, in the single scalar case, those conditions are automatically satished and do not give rise to any extra 
constraints on the functions in the held equations. It is unclear that this property persists in theories with 
multiple scalar helds, and we have not been able to determine the corresponding Lagrangian completely. 
The integrability conditions might give us some clues to accomplish this procedure. We hope to report 
the result obtained from such an approach in the near future. 
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A The ^ tensors 

In this appendix, we show the expressions for the ^ tensors introduced in Sec. 2.2 to construct Eq. (2.25). 
Expressions of and are given by Eqs. (2.22) and (2.23), and we need to construct the other ^ 
tensors appearing in Eq. (2.21), i.e., ■, and described in Sec. 2.2, 

the ^ tensors are constructed by taking the products and linear combinations of 5 “^ and that 
enjoy property S, and we need to find all such ^ tensors to construct the most general second-order tensor 
whose divergence remains of second order. In doing so, we find that not all of the ^ tensors give nontrivial 
contributions to Eq. (2.21) because there is some degeneracy and the identical terms appear from more 
than one ^ tensor. Those ^ tensors that give nontrivial contributions to Eq. (2.21) are summarized as 




abcdef _ 


+ 


I ace ^Jbdf _|_ ^Ibce^Jadf _|_ ^lade^Jbcf _|_ ^Ibde^Jacf 


^abcdef _ 

^IJ 


= CLIJKL [ € 


^aceg ^bdfh _j_ ^bceg ^adfh _|_ ^adeg ^bcfh _|_ ^bdeg ^acfh 

_|_ ^acfg^hdeh _j_ ^hcfg^adeh ^adfg^bceh _|_ ^bdfg^aceh 


Kace^Lbdf _|_ ^Kbce^Ladf _|_ ^Kade^Lbcf _|_ ^Kbde^Lacf 


biJ Qgh 


^aceg ^bdfh _|_ ^bceg^adfh _j_ ^adeg ^bcfh _|_ ^bdeg^acfh 


_j_ ^acfg^hdeh _j_ ^bcfg^adeh ^adfg^bceh _|_ ^bdfg^aceh^ 

^abcdefgh ^^^aceg^bdfh _j_ ^bceg^adfh _j_ ^adeg^bcfh _|_ ^bdeg^acfh 

_|_ ^acfg^bdeh _j_ ^bcfg^adeh ^adfg^bceh _|_ ^bdfg^aceh\ 




(A.l) 


(A.2) 


(A.3) 


where a, ajj, clijkl, b and bjj are arbitrary functions of and satisfying a/j = aji and clijkl = 
bijLK- Here we have defined as 


Jabc 


= e 


abed 1 1 

hd- 


(A.4) 


B Explicit form of Qj 

The explicit form of Qj is given as follows: 


with 


Qj = + Sf) + Sf) + + Sf) + qP + + Q^p + qP + Qp^ + Qf\ 


QP = Aj, 


(B.l) 


(B.2) 
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Qf) = -2Bij^kX^^ - Bij,KL^f, + Bij(t){,\\ (B.3) 

Qf) = (B.4) 

Qf ^ = DjKL,l5fjCj>{^cl>^\U\e^^ + 2DjjK,LX-^^cPf/ + D u K 

+ DijK,LMSb}4>\d^^^'^^^\c'^^^^ ~ B>ijK5l'f(t)\}^(l)fJ^ - -DijKSlj(l)'l^(t)^'^’'R^i (B.5) 

+ (B- 6 ) 

'^, (B.7) 

+ 2HijKL,mX'^^ 5jh(t>fe^ ^\g^ + ‘2‘HljKL,M^Jh^c^^'''^^fe^ ^\g^ 

- BljKL5lfj^(t)\^(t)fJ^(t)^''^Rgi^^ + HuKLMN^lfh^^ (B- 8 ) 
q7 = I,i6TfRj^, (B.9) 

QY^ = + 2Ju,KXJ^5)iR,J^ + 2 Jjj,k6%c^{J^\^R^/^ 

+ JlJ,KLSl%^f,^^'^%Reg^'^ ” JljKfh^J'^eg^^ (B-IO) 

Qf^ = 57 ^ 77 ?,/" - \Kj,jK5':,}trnK\%^^Rgi - lKj6^}LRce^^R,l (B-ll) 


C Euler-Lagrange Equations 

The explicit forms of the Euler-Lagrange equations from C-i-C-i are given as follows: 


7 >ab 


■pah 


(£ 1 ) = mY} 

(£ 2 ) = 2 m7j - 2 (^M 


JK 


y + 2 X'V'X’V| 


(C.l) 


,?’ + 2Af‘£,X-"'') 9'(“A?>VfJ'' 


■13 + 2Mf’„,X'“-) s''‘-4>,“3- Mg^,9‘<‘‘i,37402|.^K-|/^L|. 


- ( 

f(2) 


-lab 




(C.2) 


y - yyv y, gy t yy yj t yy yy j 

+ l^B’uc, {sH!i%<i>-’'‘‘’i>i‘e'’ + 2^'L‘<“a 2V3‘‘) - 9'<“4f]r34'*'''''<'20"'L|^i'‘ 

0,(3) „.,(3) Ax“9''“41?40"'Lfl'‘ + 2A;fT„,.Y'Ax2'f9^^^ 


-L ( 

I ^^^ILM,JK 


— 2M^^^ 

^■^'^^I[JM],KL 


(C.3) 

£“‘(£ 4 ) = -A41L2f"9”‘“A2A2<'^-R»l'“” - AA;5 x"<,'1“4,‘>“/R./a - MLi,'l'‘i2;r4'^>'''''-«./'‘ 
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(A) 


E'^’’ (Ce) 

E^^ (£ 7 ) 


- ? rAf'-*’ 


I,JK,LM^'^ + 

- 21 M>*i + Mi;:,j.LX 


“S' + X" 9'<“A,‘|'/iie/'‘ 


+ 

+ 2 
+ 2 
+ 2 


MjJ + ( 2 M, 


r(5) 

^IK,JL 


- M 


(5) 


r(5) 


+ (M, 


r(5) 


•^KL,IJ 


AM 


IJ,KL 
(5) 


X 


XL 


7X,JL 






_L9A/f^ 




-I- 


r(5) 

IM,JN,KL 


- M 


(5) 


IJ,KL,MN 


^<].J - XiTIk + 2 ( 2 AfS>,,MK- - AfSS„^) X“'l 9'<“C4'<''''''<'^ 


- - 4 (a/<5kj + A/S1./,4 - 2A2S,4i) 


M 


(7) 


+ M 


(7) 


+ M 


(7) 




^ y^^IJKLM ^ ^^^KLMIJ ^ ^^^KLIMJ) ^ 9 ^idf 9\c 9 ^\e 

^''^IJKLN,OM ^ ^^'^^IJNLK,OM ^ ‘^^^^INKLO,JMJ ^ 9 (p^^(p cp (p^g 

4 - t^LM 4J|(i iX|/ I Q/\/f(7) 

+ ^^^-‘nLIJM,KO^ ^ 9 <)zd/ <?>|c<P 9\e + ^^^hjKLM,N^ ^ 

Eo\MjjkLM,N ‘^^^hj{NK)M,L^ ^^hjMNK,L) X (p '■ <p 


\rMN ^ab 

9 


(C.4) 


(C.5) 

(C. 6 ) 


(C.7) 


D Construction of Lagrangian: single scalar-field case 

In this appendix we briefly review the construction of the Lagrangian for the most general single scalar- 
tensor theory with second order equations of motion (see Ref. [9] for more detail). The most general 
second order field equations in the single scalar-field case are given by 

Gl = A5t + {-2F'' + i + 2 D'x') (/.l>|b + (- 2 X' + 2 DX) 5Z(t>\t 

+ D4s;0|,0iV|f + \xs'^Rj> + L4s;0|f4,|f 

+ + 2-2C/1^io^'‘‘'#|f'#£ + xs'S’JiR.A + 0|>lf0l3 (D.i) 

where 0 is a scalar field and X = — {d4>)^ /2. A prime “ / ” and a dot “ ’ ” denote derivatives with respect 
to (p and X, respectively, and A,D,J, and K are arbitrary functions of (p and X, while T is related to 
the other functions as 

E = j ( 2 J -2K' + AjX^dX + W, (D.2) 

where W is an arbitrary function of (p. As explained in the main text, candidates for the Lagrangian can 
be guessed from the trace of the held equations as 

£1 = (D.3) 
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C2 = V^g , (D.4) 

£3 = (D.5) 

A = ^/^ (M<‘‘)i|»<ii,/'> + , (D.6) 

£5 = 7 ^ , (d.7) 

£6 = ^/^M( 6 ). (D. 8 ) 


By comparing the most general field equations (D.l) with the Euler-Lagrange equations obtained from 
(D.3)-(D.8), it can be seen that the free functions are related as 


A = , (D.9) 

D = - 4M(2)' + + M^^^X - 2M(^)" + + iM^^'^'X, (D.IO) 

X = - 2M(2)x - 2M(‘^)'X + 2M(5)X + 4M®X2, (D.ll) 

J = (D.12) 

K = (D.13) 

-2X" + i + 2D'X = + 2M(2)" - M(3)'X + (D.14) 

-2X' + 2DX = -M^^'>X - 2 + 2M^'^^x'^ + 3M^^'>X + 2M^^^X'^. (D.15) 

In order to identify the Lagrangian, we need to solve the above equations for First, by 

integrating Eq. (D.13) we obtain 

M(^) = -y^dX (D.16) 

Substituting Eq. (D.16) into Eqs. (D.ll) and (D.12), we then find 


= -‘^ j dX, 

= - f ^dX. 


(D.17) 

(D.18) 


Finally, we can integrate Eqs. (D.9), (D.IO), (D.14), and (D.15) to identify and with 

the help of Eqs. (D.16), (D.17) and (D.18), leading to 


and 



Af(3) = -2j ^dX, 


(D.19) 

(D.20) 


Mi^) = 2A + - 2M®'. 


(D.21) 
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